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Cosmic-rays with energies exceeding 10 19 eV are referred to as Ultra High Energy Cosmic Rays 
(UHECRs). The sources of these particles and their acceleration mechanism are unknown, and 
for many years have been the issue of much debate. The first part of this review describes the 
main constraints, that are implied by UHECR observations on the properties of candidate UHECR 
sources, the candidate sources, and the related main open questions. 

In order to address the challenges of identifying the UHECR sources and of probing the physical 
mechanisms driving them, a "multi-messenger" approach will most likely be required, combining 
electromagnetic, cosmic-ray and neutrino observations. The second part of the review is devoted 
to a discussion of high energy neutrino astronomy. It is shown that detectors, which are currently 
f"^ ■ under construction, are expected to reach the effective mass required for the detection of high energy 

£\j ' extra-Galactic neutrino sources, and may therefore play a key role in the near future in resolving 

the main open questions. The detection of high energy neutrinos from extra-Galactic sources will 
not only provide constraints on the identity and underlying physics of UHECR sources, but may 
furthermore provide information on fundamental neutrino properties. 
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I. INTRODUCTION 



Cosmic-rays (CRs) with energies exceeding ~ 10 19 eV are referred to as Ultra High Energy Cosmic Rays (UHECRs). 
The sources of these particles, which are probably extra-Galactic, and their acceleration mechanism are unknown, and 
for many years have been the issue of much debate (e.g. [3(| [H, EH III] and references therein). The first part of this 
chapter, § [H] describes the main constraints that are implied by UHECR observations on the properties of candidate 
UHECR sources. The constraints derived under the assumption that UHECRs are protons, which is supported by 
> , most observations but questioned by some (see § III Af -ij III CI § IIV|) , are summarized in § III Ft In § III Gl it is shown 
that GRBs are the only known type of sources that satisfy these constraints. Testable predictions for the spectrum 
and arrival direction distribution of UHECRs, made by the GRB model of UHECR production, are also described. 

The challenges of identifying the UHECR sources, and of p robing the physical mechanisms driving them, may be 
met with the help of high energy neutrino detectors [l|| [H, Ha ■ This is discussed in § 11111 In § 1111 Al it is shown that 
ly-j ■ detectors, which are currently under construction, are expected to reach the effective mass required for the detection 
ly-j \ of high energy extra-Galactic neutrino sources (see Figs. [9J, and may therefore play a key role in the near future 
7-H ■ in resolving the main open questions. GZK and GRB neutrinos are discussed in § IIIIBI and § MI CI respectively. In 
t-H ' § IIII D| we point out that the detection of high energy neutrinos from extra-Galactic sources will not only provide 
constraints on the identity and underlying physics of UHECR sources, but may furthermore provide information on 
f*^ ' fundamental neutrino properties. 

The main open questions associated with the production of UHECRs are summarized in § IIVI It is argued that a 
"multi-messenger" approach, combining electromagnetic, cosmic-ray and neutrino data, would be required in order 
to provide answers to these questions. 

X 

II. WHAT WE (DON'T) KNOW ABOUT THE SOURCES OF UHECRS 

The origin of CRs of all energies is still unknown (see 0, [34], [ZSj] for reviews). The cosmic ray properties change 
qualitatively as a function of particle energy, as illustrated in Fig. [1] The spectrum steepens around ~ 5 x 10 15 eV 
(the "knee") and flattens around 5 x 10 18 eV (the "ankle"). Below ~ 10 15 eV, the cosmic ray s are thought to originate 
from Galactic supernovae. However, this hypothesis has not yet been confirmed (e.g. [3 71 ] and references therein). 
The composition is dominated by protons at the lowest ener gies , and the fraction of heavy nuclei increases with 
energy. The proton fraction at ~ 10 15 eV is reduced to ~ 15% |3lT l36j. At yet higher energies, there is evidence that 
the fraction of light nuclei increases, and that the cosmic-ray flux above 5 x 10 18 eV is again dominated by protons 
[HI EH, EE] • The composition change and the flattening of the spectrum around 10 19 eV suggest that the flux above 
and below this energy is dominated by different sources. At energies of £19 = £7/10 eV ~ 1 the Larmor radius of 
CRs in the Galactic magnetic field is 

Rh = l^B^ ZBllsEwZ-^pc, (1) 
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FIG. 1: A schematic description of the differential CR spectrum, dJ/dE, with some comments on what we know (or don't) 
about the composition and origin of the CRs. 



where B = 10 .B_5.5/zG is the value of the Galactic magnetic field and Z is the CR charge. Since the Galactic 
magnetic field cannot confine protons above 10 19 eV, it is believed that the nearly isotropic cosmic ray flux at 
E > 5 x 10 18 eV originates from extra-Galactic (XG) sources. In what follows we focus on this XG component. 



A. Composition 

At low energy, < 10 TeV, CR particles are detected by space or balloon born detectors, which provide a direct 
measurement of the primary CR composition. At higher energies, the flux is too low to be detectable by space/balloon 
born detectors, and CRs are detected indirectly through the "air-showers" , the large number of lower energy particles, 
they produce as they propagate and lose energy in the atmosphere. The low flux at the highest energies, 

J(> 10 20 eV) w 1/100 km 2 yr27rsr, (2) 

requires detectors with effective area of many 100's km 2 . The primary composition is constrained at high energies 
mainly by the average and variance of X max , the depth in the atmosphere at which the shower contains the largest 
number of high energy particles, obtained for showers of fixed energy (fluctuations in individual shower development 
are large, leading to fluctuations in the depth of maximum which are not small compared to the dependence on the 
primary mass). X max is larger for higher energy particles. Since a high energy heavy nucleus behaves roughly as a 
group of independent lower energy nucleons, X max and its variance are larger at fixed energy for lighter nuclei. 

Fig. [2] presents the main evidence for the transition to lighter nuclei at higher energy: X max grows with energy 
faster than model predictions for fixed composition, becoming consistent with pure proton composition at ~ 10 18 eV. 
At the highest observed energies, there is some discrepancy between the results reported by the HiRes observatory 
and by the Pierre Auger Observatory (PAO). While the HiRes observatory reports the average and variance of X max 
to be consistent with a pure proton composition all the way up to 10 197 eV (Fig. [3]), the PAO reports X max and ax 
evolution which suggests a transition back to heavier nuclei at the highest energies (Fig. SJ. 

Th e origin of this discrepancy is not yet understood. However, a few comments are in place. It was noted in 
\l0l\ that the analysis of the PAO data, presented in Q, is not self consistent: according to this analysis, ax 
measured at the highest energy implies an Fe fraction > 90%, while the measured value of < A" max > implies an Fe 
fraction < 60%. This inconsistency may reflect some experimental problem, but may also reflect a modification of 
the hadronic interaction cross section which is not accounted for in the models used for shower calculations. It should 
be emphasized that the theoretical AT max calculations depend on extrapolation of hadronic models to energies well 
beyond those currently tested in accelerators. The theoretical and experimental uncertainties in the extrapolation of 
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FIG. 2: Average depth of shower maximum as function of energy: Measurements by the HiRes detector compared to predictions 
for proton and iron primaries based on various model extrapolations of the pp cross section. Adapted from [3| . 
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FIG. 4: Same as Fig. [5] for the PAO data. Adapted from [f|. 
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FIG. 5: Left: A comparison of direct numerical calculations of the effective CR life time (solid lines) with the analytic 
approximation of HJ using {E c , ep = 9.1 x 10 18 eV, r , ep = 0.5 x 10 9 yr, E c , n = 3.5 x 10 20 eV, To,* = 1.4 x 10 7 yr} (dashed 
line), for CR generation following dh/dE(E, z) oc (1 + z) m E~ a . Right: The local (z = 0) energy generation rate as measured 
by Auger [35l | and Hires [I| assuming that the CRs are purely protons, for a — 1 = 1 (For different values of a, the spectrum 
should be multiplied by an energy independent factor (a — 1); Q = h). Statistical and systematic errors in the experimental 
determination of event energies lead to ~ 50% uncertainty in the flux at the highest energies. The absolute energy scales of 
the Auger and Hires data where not altered in this figure. Adapted from [60|. 



the pp cross-section to center-of-mass energies > 100 TeV are a possible source of biases in shower reconstruction (e.g. 
89]). It is therefore difficult to draw a firm conclusion regarding primary composition at the highest energies based 
on current shower measurements. 



B. Generation rate & spectrum 

Let us assume first that the UHECRs are protons of extra-Galactic origin. As they propagate, high-energy protons 
lose energy as a result of the cosmological redshift and as a result of production of pions and e + e— pairs in interactions 
with cosmic microwave background (CMB) photons. The local intensity of UHECRs may be written as 

where dho(E) / dE is the local (z = 0) proton production rate (per unit volume and proton energy) and t s. is the 
effective energy loss time of the proton (this equation is, in fact, a definition of t c ff.). The left panel of Fig. [3] shows 
t e s. for proton generation following dn/dE(E, z) oc (1 + z) rn E~ a . The rapid decrease in the effective life time, or 
propaga tion distance ct c g. , above ~ 6 x 10 19 eV, commonly termed the "Greisen-Zatsepin-Kuzmin (GZK) suppression" 
[471 Il09j | , is due to photo-production of pions by the interaction of protons with CMB photons (The proton threshold 
energy for pion production on ~ 10 -3 eV CMB photons is ~ 10 20 eV). Since proton propagation is limited at high 
energies to distances <C c/Hq, e.g. to ~ 100 Mpc at 10 20 eV, the dependence of t c g. on redshift evolution (m) is not 
strong. 

Using Eq. © and the measured UHECR intensity, it is straightforward to infer the local production rate of UHECRs. 
The right panel of figure [5] shows that the energy generation rate above 10 19 ' 5 eV is roughly constant per logarithmic 
CR energy interval, a ~ 2 and 

E2 dh^E)_ ^ 10 43.5 erg/Mpc 3 yr (4) 

dE 

In other words, the observed CR spectrum is consistent with a generation spectrum dh/dE oc E~ 2 modified by the 
GZK suppression. Since both observations and models for particle acceleration in collisionless shocks, which are 
believed to be the main sources of high energy particles in many astrophysical systems, typically imply a ~ 2 (see 
[H], [98j for reviews of particle acceleration in non-relativistic and relativistic shocks respectively) , this supports the 
validity of the assumption that UHECRs are protons produced by extra-Galactic objects. 

The following point should, however, be made here. Heavy nuclei lose energy by interaction with CMB and IR 
photons, that leads to spallation. Since the effective life time of such nuclei is not very different from that of protons, 



5 



the consistency of the observed spectrum with a model of extra-Galactic sources of protons with generation spectrum 
of dh/dE oc E~ 2 can not be considered as a conclusive evidence for the UHECRs being protons. 

One of the important open questions is at what energy the transition from Galactic to extra- Galactic (XG) sources 
takes place. A simple model with E 2 dho(E)/dE = 5 x 10 43 erg/Mpc 3 yr and a transition from Galactic to XG sources 
at 10 19 eV is consistent with observations [6CJ. In such a model, the Galactic flux is comparable to the XG one 
at 10 19 eV, and negligible at > 10 19 ' 5 eV. Other models, however, have been proposed, in which the Galactic- XG 
transition occurs well below 10 19 eV (e.g. [30| and references therein). Such models are motivated mainly by the 
argument that they allow one to explain the ~5x 10 18 eV spectral feature by pair production (in proton interactions 
with the CMB). The transition energy in such models is therefore well below 5 x 10 18 eV. As explained in [6(J, a 
Galactic-XG transition at ~ 10 18 eV requires fine tuning of the Galactic and XG contributions (to produce the smooth 
power- law observed) , and is disfavored by the data: it requires that Auger systematically underestimates the energy 
of the events by 40% (well above the stated uncertainty) and it requires dhp^c/de oc e~ 2 - 7 , which is inconsistent with 
the > 10 19 eV data. 

Finally, one notes that if the generation spectrum of XG CRs extends over many decades below 10 19 eV, the 
total XG CR energy production rate, Qxq, might exceed significantly the UHECR production rate, Q^il c y = 
(E 2 dh /dE) 

>io 19 eV; given by Eq. Q. For the dh/dE oc l/E 2 spect rum inferred from observations, the "bolometric 
correction" will be Qxg/Qio 19 oV ~ ln(10 20 eV/£ ; m i n ) ~ 10, where -Emm <C 10 19 eV is the low energy to which the 
spectrum extends. 



C. Anisotropy: Source and composition clues 



The propagation of UHECRs is limited at the highest energies to distances ~ 100 Mpc. The galaxy distribution 
is not homogeneous over such a distance scale. Thus, if the distribution of UHECR sources is correlated with that 
of galaxies, one expects an anisotropy in the UHECR arrival direction distribution reflecting the inhomogeneity of 
the galaxy distribution [l03j |. Fig. [5] shows the integrated galaxy density out to 75 Mpc and the predicted anisotropy 
of the UHECR intensity. Also shown are the (angular) positions of the 27 Auger events with energy exceeding 
5.7 x 10 19 eV. The distribution of these events is inconsistent with isotropy at a 98% confidence level for a source 
density n s = 10~ 4 Mpc~ 3 , corresponding to the lowest allowed source density (see § III D|) . and at a 99% confidence 
level for n s = 10 _2 Mpc~ 3 , corresponding to the density of galaxies. The angular distribution of CR arrival directions 
is consistent with a UHECR source distribution that follows the galaxy distribution (for detailed discussions see 
88]). This provides some support to the association of the sources with known extra-Galactic astrophysical 
objects llo| . The more recent PAO analysis (valid for n s — > oo) of a larger number of events, 58 above 5.5 x 10 19 eV, 



yields inconsistency with isotropy at a 99% confidence level 

UHECRs may suffer significant deflections as they cross dense large scale structures, such as galaxy clusters and 
large scale galaxy filaments, in which the energy density of the plasma is large enough to support strong magnetic 
fields. Such deflections may distort the anisotropy pattern expected based on the galaxy distribution. An estimate of 
the expected deflection may be obtained assuming that all large scale structures support a magnetic field with energy 
density comprising a fraction eg of the plasma thermal energy density. The deflection expected in this case for a 
propagation distance d is (see 59, 62] for a detailed derivation) 



1 Mpc V0.1 100 Mpc 10 kpc/ V0.01/ V lo2 ° eV 

Here, Z is the particle charge, / is the fraction of the volume filled by filaments of diameter L, and A is the field 
coherence length. The deflections are not expected therefore to distort significantly the anisotropy map. 

The anisotropy signal provides also a test of the primary UHECR composition. If one records an anisotropy signal 
produced by heavy nuclei of charge Z above an energy -Ethr, one should record an even stronger (possibly much 
stronger) anisotropy at energies > E t hr/Z due to the proton component that is expected to be associated with the 
sources of the heavy nuclei. This is due to the fact that particles of similar rigidity E/Z propagate in a similar manner 
in the inter-galactic magnetic field and based on the plausible assumptions that (i) a source accelerating particles of 
charge Z to energy E will accelerate protons to energy E/Z, and (ii) there are at least as many protons accelerated as 
there are heavy nuclei. The anisotropy signal is expected to be stronger at lower energy since the signal increases as the 
number of particles produced by the source, E~ a+1 , while the background increases as the square- root of the number 
of all observed CRs, ~ £ ,_ ( 2 - 7_1 )/ 2 (see [69| for a detailed discussion). Thus, if the PAO > 5.7 x 10 19 eV anisotropy 
signal is real, the lack of detection of stronger anisotropy at lower energy disfavors a heavy nuclei composition at 
- 6 x 10 19 eV. 



FIG. 6: Top: The integrated galaxy density out to a distance of 75 Mpc, normalized to the mean integrated density. The 
contours are logarithmic, ranging from 0.5 to 4 with three contours per density doubling. Dashed curves represent under- 
density. Bottom: The positions of the 27 Auger events with energy exceeding 5.7 x 10 19 eV Q, overlaid on the UHECR 
intensity map, J(fi), predicted in a model in which the UHECR source distribution follows the galaxy density distribution 
(with a bias b[S] = 1 + 5 for S > 0, b = otherwise, where S is the fractional galaxy over density). The coordinates are Galactic 
and J is normalized to its all sky average. The contours denote J/J = (0.7,0.9, 1, 1.1, 1.3, 1.5), with dashed lines representing 
under-density. The thick solid line denotes the super-galactic plane. The dashed-dotted green line marks the boundary of 
Auger's coverage (corresponding to a zenith angle of 60°). Adapted from [591 ]. 



The arrival directions of the 27 PAO events and ~ 30 HiRes events above 6 x 10 19 eV show no evidence for 
"repeaters" , i.e. multiple events that may be associated with a single source given the small deflection angles expected. 
The lack of repeaters implies that the number of sources contributing to the flux, N s , should satisfy N s > N 2 , where N 
is the number of events (for identical sources each producing on average N/N s events and N 2 /N s <C 1, the probability 
for repeaters is ~ N 2 /N s ). This suggests that there should be more than ~ 10 3 5 independent sources contributing 
to the (all sky) flux (note that HiReS and PAO observed the northern and southern hemispheres respectively). For 
protons, the effective propagation distance is ~ 200 Mpc, see Fig. [5J implying a lower limit on the source density of 



D. Source density 



n s > 10~ 4 Mpc 



-3 



(6) 



(for a more detailed analysis see [39l. fl03]). For comparison, the density of galaxies is roughly 10 



Mpc 
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E. Source constraints: Minimum power and speed 



The essence of the challenge of accelerating particles to > 10 19 eV can be understood using the following simple 
arguments ([97|], for a more detailed derivation see 93]). Consider an astrophysical source driving a flow of magnetized 
plasma, with characteristic magnetic field strength B and velocity v. Imagine now a conducting wire encircling the 




FIG. 7: Potential drop generated by an unsteady outflow of magnetized plasma. 



source at radius R, as illustrated in Fig. [7J The potential generated by the moving plasma is given by the time 
derivative of the magnetic flux $ and is therefore given by V ~ /3BR where /3 = v/c. A proton which is allowed to 
be accelerated by this potential drop would reach energy E ~ /3eBR. The situation is somewhat more complicated 
in the case of a relativistic outflow, with T = (1 — /3 2 ) -1 / 2 ^s> 1. In this case, the proton is allowed to be accelerated 
only over a fraction of the radius R, comparable to R/Y. To see this, one must realize that as the plasma expands, its 
magnetic field decreases, so the time available for acceleration corresponds to the time of expansion from R to, say, 
2R. In the observer frame this time is R/c, while in the plasma rest frame it is R/Yc. Thus, a proton moving with 
the magnetized plasma can be accelerated over a transverse distance ~ R/T. This sets a lower limit to the product 
of the magnetic field and source size, which is required to allow acceleration to E, BR > YE/e(3. This constraint also 
sets a lower limit to the rate L at which energy should be generated by the source. The magnetic field carries with it 
an energy density B 2 /8tt, and the flow therefore carries with it an energy flux > vB 2 /8ir (some energy is carried also 
as plasma kinetic energy), which implies L > vR 2 B 2 and therefore 



Another constraint on the source results from the requirement that the acceleration is not suppressed by synchrotron 
emission of the accelerated particle. Let us consider a relativistic source. Denoting by B' the magnetic field in the 
plasma rest frame, the acceleration time of a proton is t' acc > E' /eB'c where E' = E/Y is the proton energy in the 
plasma frame. The synchrotron loss time, on the other hand, is given by t' syn » (m p /m e ) 2 (6TrE' /a T cy 2 B' 2 ) where 
7' = E'/nipC 2 . Requiring t' &cc < t' syn sets an upper limit on B' (which depends on E and Y). Requiring this upper 
limit to be larger than the lower limit derived in the previous paragraph, B' R > E/e, sets a lower limit to Y (which 
depends on R and E) . Relating the source radius R to an observed variability time (of the radiation emitted by the 
source) through R = 2Y 2 cSt, the lower limit is [93[ 

/ E \3/4 / r. s -1/4 

r>1 ° 2 UM (is;;) ■ W 

This implies that the sources must be relativistic, unless their characteristic variability time exceeds ~ 10 6 s. 



F. Summary of source constraints 



The evidence for a transition to a light composition, consistent with protons, at few xlO 18 eV (§ 111 Aj) . the consis- 
tency of the spectrum above ~ 10 19 eV with a dh/dE oc 1/E 2 generation spectrum modified by the GZK suppression 
(§ HI Bp . and the hints for a light composition from the anisotropy signal (§ III Cp . suggest that the UHECRs are 
protons produced by extra-Galactic sources. If this is indeed the case, the discussion of the preceding sections implies 
that their sources must satisfy several constraints: 

• The sources should produce protons with a local (z = 0) rate and spectrum (averaged over space and time) 
E 2 dh /dE 10 43 - 5 erg/Mpc 3 yr; 

• The density of sources (contributing to the flux at ~ 5 x 10 19 eV) should satisfy n s > 10 _4 Mpc -3 ; 

• The power output of the individual sources should satisfy L > 10 45 ' 5 r 2 /3 _1 erg/s; 

• The Lorentz factor of the flow driven by the source must satisfy Y > 10 2 ((5t/10ms) -1 / 4 where St is the charac- 
teristic source variability time. 
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No sources that satisfy the constraint L > 10 46 erg/s are known to lie within a ~ 100 Mpc distance. One may 
argue, of course, that there are "dark sources" , i.e. sources that produce such power output (and UHECRs) but do 
not produce much radiation and are hence not known. One can not rule out the existence of such sources. On the 
other hand, we do not have direct evidence for their existence either. Putting aside such a caveat, the lack of known 
sources of sufficient luminosity suggests that the sources are transient. The transient duration T must be shorter than 
the time delay between the arrival of photons and protons from the source. The protons are delayed due to magnetic 
field deflection by At ~ 8 2 d/c, where 9 is estimated in Eq. ([S]). This yields 

A ^ d )~ lo4 (lM^) 2 (io4v)" 2yr - (9) 

Due to the random energy loss of the protons during their propagation, and due to the possibility of multiple paths 
between source and observer, the arrival of protons of energy E is delayed and spread over a similar time At(E, d). 
For T < At, the effective number density of sources contributing to the flux at energy E is ~ n s At[E, d e g(E)], where 
h s is the transient rate (per unit volume) and d c fi(E) ~ ci e ff (-£■). 



G. "Suspects", Predictions 

Only two types of sources are known to satisfy the above minimum power requirement: active galactic nuclei ( AGN ) 
- the brightest known steady sources, and gamma-ray bursts (GRBs) - the brightest known transient sources [ 1 1 1| . 
Both AGN (e.g. [83() and GRBs [73|, |90, |93( have therefore been suggested to be UHECR sources. The absence 
of AGN with L > 10 46 erg s _1 within the GZK horizon had motivated the suggestion j44j that UHECRs may be 
produced by a new, yet undetected, class of short duration AGN flares resulting from the tidal disruption of stars or 
accretion disk instabilities. The existence of tidal disruption flares is likely. However, they are yet to be dete cted and 
whether their properties are consistent with the constraints derived above is yet to be determined (see also jl04j ). 

Let us consider then the GRB transients. First, consider the minimum power and minimum speed constraints 
that should be satisfied by individual sources: Eq s. Q and ©. For GRBs, the (luminosity function averaged) peak 
luminosity is L 1 w 10 52 erg/s (HO, [§2, note that (50| gives L R n-3no k R y which is « 0.1 of Lq.i-io McV given in 92]), 
and typical values of T and St are V ~ 10 25 and St ~ 10 ms [l0, UM, US [§3 ■ Thus, both constraints are satisfied. It is 
worth noting that T > 10 2 is inferred for GRBs based on the photon spectrum (in order to avoid large pair production 
optical depth), i.e. based on arguments which are different than those leading to the T > 10 2 constraint of Eq. (|SJ). 

Next, let us consider the global constraints on the rate, Eq. ([6]), and average energy production rate, Eq. (j4j, of 
the sources. The local, z = 0, GRB rate is h z=0 ~ 10 _9 Mpc~ 3 yr _1 (assuming h s evolves rapidly with redshift, 
following the star formation rate, i.e n z=0 <C n z s =1 - 5 [lo|, H), implying, using Eq. ©, n s (E) ~ h z s =Q At[E, d e g (E)] — 
10- 4 (4ff/200Mpc) 2 (£;/0.5 x 10 20 eV)- 2 Mpc- 3 , consistent with Eq. ©. The local, z = 0, GRB energy production 
rate in ~ 1 MeV photons is given by n z=a L y At, where At is the effective duration (the average ratio of the fluence 
to the peak luminosity) corrected for redshift (the observed duration is 1 + z larger than the duration at the source), 
At w 10s/(l + z) ~ 4 s (using z — 1.5 as a characteristic redshift). This yields E~ = L 1 At w 10 52 5 erg and 
QmcV grb = n^Ej 10 43 5 erg/Mpc 3 yr, similar to the required UHECR energy production rate given in Eq. ([J}, 
QIqw cV = (E 2 dn /dE) >w i9 cV ks 10 43 5 erg/Mpc 3 yr (for a more detailed discussion see [13, l9rj l99l : for additional 



energy production by "low-luminosity GRBs" and "heavy baryon loading GRBs" see [76| and [106l | respectively, and 
references therein). 

As noted at the end of § III Bt if the generation spectrum of XG CRs extends over many decades below 10 19 eV, the 
total XG CR energy production rate, Qxg°' ma y exceed significantly the UHECR production rate, <5xg/Qio 19 cV ~ 
10. Estimating the ratio of Q X g° to the total photon energy production by GRBs, Q^qrb^ as Qxg /Q^GKB = 
Qxg°/'9mcV grb ~ t0 is, however, quite uncertain. This is due to uncertainties in the redshift evolution of the 
GRB rate and luminosity function, in the "bolometric correction" for the CR production rate, and in the bolometric 
correction, Qj = q-rb/QmcV grb > 1' that should also be applied to the photons. 

If GRBs are the sources of UHECRs, then some int eres ting predictions can be made regarding the spectrum and 
angular distribution of events at the highest energies [l05j ] . Due to the rapid decrease of d e s with energy, the total 
number of sources contributing to the flux, ~ (4Tr/3)dl e n s At[E, d c g{E)] drops rapidly with energy. This implies that, 
for n s ~ 10 _9 Mpc~ 3 yr~ 1 and adopting the estimate of Eq. (J5J) for the deflection angle, only a few sources contribute 
to the flux above ~ 3 x 10 20 eV. Moreover, the spectrum of these sources should be rather narrow, AE/E ~ 1, since 
the energy dependent time delay At(E, d) implies that higher (lower) energy particles arrived (will arrive) in the past 
(future). Testing this prediction, which requires a large number of events detected above ~ 3 x 10 20 eV, may require 
large exposure, exceeding even that of PAO, which may be provided by space born detectors [8ll[87|. 
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FIG . 8: The upper bound imposed by UHECR observations on the extra-Galactic high energy muon neutrino + intensity 
[24l . Il0l| (lower-curve: no evolution of the energy production rate, upper curve: assuming evolution following star formation 
rate), compared with the atmospheric muon-neutrino background and with several experimental upper bounds (solid lines). The 
theoretical bound does not include the effect of neutrino oscillations. Such oscillations are expected to change the i/ e : : v T 
flavor ratio from f : 2 : to f : 1 : 1 (e.g. [H3]), leading to an upper bound which is ~ 1/2 that shown in the figure for each flavor. 
Shown are the muon and all flavor upper bounds of the optical Cerenkov observatories AMANDA [ToL [TO] and BAIKAL [22| . 
the all flavor upper bounds of the coherent Cerenkov radio detectors RICE [63|] and ANITA [4^], and the v T upper bound of the 
PAO [|. The curve labelled "GZK" shows the muon neutrino intensity (not corrected for oscillations) expected from UHECR 
proton interactions with micro- wave background photons [2S|. Black dashed curves show the expected sensitivity (for few 
years operation) of 0.1 Gton (ANTARES, http://antares.in2p3.fr/) and 1 Gton (IceCube, http://icecube.wisc.edu/; Km3Net, 
http://www.km3net.org/home.php) optical Cerenkov detectors. The blue dashed curve is the expected sensitivit y o f detectors 
of few 100 Gton (few 100 km 3 ) effective mass (volume), that may be achieved with proposed radio detectors [l^. l26l. [271. [66p or 
with proposed (optical) extensions of IceCube [HiJ. For a detailed discussion of the current experimental status see fig. |56|]. 

III. HIGH ENERGY NEUTRINO ASTRONOMY 

UHECR sources are likely to be sources of high energy neutrinos. The interaction of high energy protons (nucleons) 
with radiation or gas, either at or far from the source, leads to production of charged pions, via pj and pp(n) 
interactions, which decay to produce neutrinos (e.g. p + 7 — > n + 7r + , tt + — > fi + + — > e + + + + u e ). In § IIII Al 
we estimate the minimum detector size, which is required to detect such neutrinos. In § IIII Bl we comment on the 
importance of the detection of "GZK neutrinos". The prospects for detection of GRB neutrinos, and the possible 
implications of such detection for the study of GRBs, are discussed in § 1111 CI Prospects for the study of fundamental 
neutrino properties using high energy GRB neutrinos are discussed in § IIII Dl For most of the discussion of this 
section, we adopt the assumption that UHECRs are protons. 

A. Neutrino flux upper bound, Detector size, Detectors' status 

The energy production rate, Eq. Q, sets an upper bound to the neutrino intensity produced by sources which, as 
GRBs and AGN jets, are for high-energy nucleons optically thin to pj and pp(n) interactions. For sources of this 
type, the energy generation rate of neutrinos can not exceed the energy generation rate implied by assuming that all 
the energy injected as high-energy protons is converted to pions (via pj and pp{n) interactions). Using Eq. the 
resulting upper bound (y^ + D^, neglecting mixing) is pil Il0l| 

^ < - ^ GoCCJ G-c-V SI -.. (1 0, 

Here tn is the Hubble time and the 1/4 factor is due to the fact that charged and neutral pions (which decay to 
photons) are produced with similar probability, and that muon neutrinos carry roughly half the energy of the decaying 
pion. In the derivation of Eq. (|10|) we have neglected the redshift energy loss of neutrinos produced at cosmic time 
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t < tii, and implicitly assumed that the cosmic-ray generation rate per unit (comoving) volume is independent of 
cosmic time. The quantity £z in Eq. [TU] has been introduced to describe corrections due to redshift evolution and 
energy loss. 

The upper bound is compared in Fig. [5] with the current experimental limits and with the expected sensitivity of 
planned neutrino telescopes. The figure indicates that km-scale (i.e. giga-ton-scale) neutrino telescopes are needed 
to detect the expected extra-Galactic flux in the energy range of ~ 1 TeV to ~ 1 PeV, and that much larger effective 
volume is required to detect the flux at higher energy. The Baikal, AMANDA, and ANTARES optical Cerenkov 
telescopes have proven that the construction of km-scale neutrino detectors is feasible, and the IceCube detector, the 
construction of which is well underway, is expected to reach its designed target effective mass of ^ 1 Gton in 2011. 

B. GZK neutrinos 

As discussed in § III Bt protons of energy exceeding the threshold for pion production in interaction with CMB 
photon, ~ 5 x 10 19 eV, lose most of their energy over a time short compared to the age of the universe. If UHECRs 
are indeed protons of extra-Galactic origin, their energy loss should produce a neutrino intensity similar to the upper 
bound given by Eq. (|T0|) . Since most of the pions are produced in interactions with photons of energy corresponding 
to the A-resonance, each of the resulting neutrinos carry approximately 5% of the proton energy. The neutrino 
background is therefore close to the bound above ~5x 10 18 eV, where neutrinos are produced by ~ 10 20 eV protons. 
The intensity at lower energies is lower, since protons of lower energy do not lose all their energy over the age of the 
universe (The GZK intensity in figure [8] decreases at the highest energies since it was assumed that the maximum 
energy of protons produced by UHECR sources is 10 21 eV). The results of detailed calculations of the expected GZK 
neutrino intensity 143:] are in agreement with the qualitative analysis presented above. 

The detection of GZK neutrinos will be a milestone in neutrino astronomy. Most important, it will allow one to test 
the hypothesis that the UHECRs are protons (possibly somewhat heavier nuclei) of extra-Galactic origin (e.g. [75| 
and references therein). Moreover, measurements of the flux and spectrum would constrain the redshift evolution of 
the sources. Finally, detection of ultra-high energy neutrinos may allow one to test for modifications of the neutrino 
interaction cross section due to new physics effects at high (100 TeV) energies [l5l l27l. l65j. 

C. Neutrinos from GRBs 

GRB gamma-rays are believed to be produced within a relativistic expanding wind, a so called "fireball" , driven 
by rapid mass accretion onto a newly formed stellar-mass black hole. It is commonly assumed that electrons are 
accelerated to high energy in collisionless shocks taking place within the expanding wind, and that synchrotron 
emission from these shock accelerated electrons produces the observed 7-rays (see [70, 0, H] for reviews). If 
protons are present in the wind, as assumed in the fireball model, they would also be accelerated to high energy in 
the region were electrons are accelerated. If protons are indeed accelerated, then high energy neutrino emission is also 
expected. 

1. 100 TeV fireball neutrinos 

Protons accelerated in the region where MeV gamma-rays are produced will interact with these photons to produce 
pions provided that their energy exceeds the threshold for pion production, 

E 1 E w 0.2r 2 GeV 2 . (11) 

Here, E-y is the observed photon energy. The T 2 factor appears since the protons and photon energies in the plasma 
rest frame (where the particle distributions are roughly isotropic) are smaller than the observed energy by the Lorentz 
factor T of the outflow. For T « 10 2 5 and E y = 1 MeV, proton energies ~ 10 eV are required to produce pions. 
Since neutrinos produced by pion decay typically carry 5% of the proton energy, production of ~ 10 eV neutrinos 
is expected |100| . 

The fraction of energy lost by protons to pions, f v , is f n ~ 0.2 jHlL Il00| . Assuming that GRBs gen erate the observed 
UHECRs, the expected GRB muon and anti-muon neutrino flux may be estimated using Eq. (|T0|) [lOfjj Il0l| . 

« 10- 8 A ( -J^-gg; ) GeVcm-VV- 1 . (12) 
v 0.2 V 10 44 erg/Mpc 3 yry v ' 
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FIG . 9: AMANDA flux upper limits (solid lines, 90% confidence) for muon neutrino energy spectra predicted by the models of 
[73,JloO] for the ~ 100 TeV internal shock fireball neutrino s (§IIIIC1[) , and for the muon neutrino energy spectrum predicted 
by |85l] for the prec ursor supernova ( "supranova" ) model (S IIII C 2b - The upper bounds are compared with the fluxes predicted 
by the models f [95l . ll00l ]- thick dotted line, thin dotted line. [85|- dot-dashed line). Adapted from 



This neutri no spectr um extends to ~ 10 16 eV, and is suppressed at higher energy due to energy loss of pions and 
muons f&il llOOl llOlj (for the contribution of Kaon decay at high energy see pj|). Eq. (fl"2")l implies a detectio n rat e 
of ~ 10 neutrino-induced muon events per year (over 4tt sr) in a 1 Gton (1 cubic-km) detector (13, I2H |49|. [77l . Il00j |. 
The upper limit on the GRB neutrino emission provided by the AMANDA (~ 0.05 Gton) detector approaches the 
flux predicted by Eq. (|T2"j) . see Fig. [51 and the 1 Gton IceCube detector, which will be completed at the beginning of 
2011, will reach a sensitivity that may allow one to test this model's predictions [l|. 

Since GRB neutrino events are correlated both in time and in direction with gamma-rays, their detection is practi- 
cally background free. The main background is due to atmospheric neutrinos, which produce neutrino-induced muons, 
travelling in a direction lying within a cone of opening angle A9 around some direction, at a rate 



4 x 10 
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AO 

05° 



—) 

lOOTeV/ 



-0 



km 2 yr , 



(13) 



with j3 = 1.7 for E < 100 TeV and j3 — 2.5 for E > 100 TeV. At high energies, the neutrino induced muon propagates 
at nearly the same direction as the incoming neutrino, and km-scale neutrino telescopes will be able to determine the 
incoming neutrino direction to better than ~ 0.5°. For a known source direction, therefore, the neutrino search is 
practically background free. 



2. Te V neutrinos 



The 100 TeV neutrinos discussed in the previous sub-section are produced in the same region where GRB 7-rays are 
produced and should therefore accompany the 10 to 100 s 7-ray emission phase (note, however, that it was pointed out 
in [lH that if the late, ~ 10 4 s, X-ray/UV flares are produced by late internal shocks within the fireball, the emission 
of 100 TeV neutrinos may be extended to accompany these flares). Their production is a generic prediction of the 
fireball model: it is a direct consequence of the assumptions that energy is carried from the underlying engine as kinetic 
energy of protons and that 7-rays are produced by synchrotron emission of shock accelerated particles. Neutrinos 
may be produced also in other stages of fireball evolution, at energies different than 100 TeV. The production of 
these neutrinos is dependent on additional model assumptions. We discuss below some examples of ~ 1 TeV neutrino 
emission predictions, that depend on the propertie s of t he GRB progenitor. For a discussion of ~ 10 18 eV neutrino 
emission during the afterglow phase see [zH HH Il02j | and the reviews [z3, [Zll HH ■ 

The most widely discussed progenitor scenarios for long-duration GRBs involve core collapse of massive stars. In 
these "collapsar" models, a relativistic jet breaks through the stellar envelope to produce a GRB. For extended or 
slowly rotating stars, the jet may be unable to break through the envelope. Both penetrating (GRB producing) and 



12 



"choked" jets can produce a burst of ~ 10 TeV neutrinos by interaction of accelerated protons with jet photons, while 
the jet propagates in the envelope [13, [Z2, 86] (it was pointed out in [TtJ that neutrino production by kaon decay may 
dominate over the pion decay contribution, extending the neutrino spectrum to ~ 20 TeV) . The estimated event rates 
may exceed ~ 10 2 events per yr in a km-scale detector, depending on the ratio of non-visible to visible fireballs. A 
clear detection of non- visible GRBs with neutrinos may be difficult due to the low energy resolution for muon-neutrino 
events, unless the associated supernova photons are detected. 

In the two-step "supranova" model, interaction of the GRB blast wave with the supernova shell can lead to detectable 
neutrino emission, either through nuclear collisions with the dense supernova shell or through interaction with the 
intense supernova and backscattered radiation field [42|, EH 85] . As indicated by Fig. [9] the upper limits provided by 
AMANADA on the muon neutrino flux suggest that "supranova" s do not accompany most GRBs. 

D. Neutrino physics prospects 

In addition to testing the GRB model for UHECR production and to providing a new handle on the physics of 
GRB sources, detection of high energy GRB neutrinos may provide information on fundamental neutrino properties 

[iooj . 

Detection of neutrinos from GRBs could be used to test the simultaneity of neutrino and photon arrival to an 
accuracy of ~ 1 s. It is important to emphasize here that since the background level of neutrino telescopes is very 
low, see Eq. (|13p. the detection of a single neutrino from the direction of a GRB on a time sale of months after the 
burst would imply an association of the neutrino with the burst and will therefore establish a time of flight delay 
measurement. Such a measuremen t wil l allow one to test for violations of Lorentz invariance (as expected due to 
quantum gravity effects) [H], |H, H3, [t 00]), and to test the weak equivalence principle, according to which photons and 
neutrinos should suffer the same time delay as they pass through a gravitational potential. With 1 s accuracy, a burst 
at 1 Gpc would reveal a fractional difference in (photon and neutrino) speed of 10~ 17 , and a fractional difference in 
gravitational time delay of order 10~ 6 (considering the Galactic potential alone). Previous applications of these ideas 
to supernova 1987A (see [25[ for review), yielded much weaker upper limits: of order 10~ 8 and 10 -2 respectively. 
Note that at the high neutrino energies under discussion deviations of the propagation speed from that of light due to 
the finite mass of the neutrino lead to negligible time delay even from propagation over cosmological distances (less 
than - 10~ 10 s at 100 TeV). 

High energy neutrinos are expected to be produced in GRBs by the decay of charged pions, which lead to the 
production of neutrinos with flavor ratio $„ c : : <i>„ T = 1:2:0 (here <& Vl stands for the combined flux of v\ 
and v{). Neutrino oscillations then lead to an observed flux ratio on Earth of $„ c : : $„ 7 =1:1:1 [2l|, |68( 
(see, however [581]). Up-going r's, rather than fi's, would be a distinctive signature of such oscillations. It has 
furthermore been been pointed out tha t fla vor measurements of astrophysical neutrinos may help determining the 
mixing parameters and mass hierarchy [J08] , and may possibly enable one to probe new physics (20 . |68| . 

IV. OUTLOOK: OPEN QUESTIONS & MULTI-MESSENGER ASTRONOMY 

The validity of the constraints imposed on the properties of candidate UHECR sources, as summarized in § III F[ 
depends on the validity of the inference that the highest energy particles are protons, and on the validity of the 
assumption that the particles are accelerated by some electromagnetic process, for which the constraints derived in 
§ III El are valid. The inference that the highest energy particles are protons is supported by the HiRes and PAO 
UHECR spectrum, by the properties of air showers as measured by HiRes, and by the anisotropy hints. However, the 
shower properties reported by PAO appear to be inconsistent with a pure proton composition at the highest energy 
(and possibly also with a heavy nuclei composition, see § III A[) . Given this, and the fact that the pp cross section 
at the high energies under discussion is not well known, the possibility that the highest energy particles are heavy 
nuclei can not yet be excluded. If the particles are indeed heavy nuclei of charge Z , the minimum power requirement, 
Eq. ([7]), would be reduced by a factor Z 2 , and could possibly be satisfied by local steady sources like AGN (e.g. [Hj). 

Thus, although we have strong arguments suggesting that UHECR sources are protons produced by transient XG 
sources, and that the sources should satisfy the constraints given in § III Ft which point towards GRBs being the 
likely sources, we are still missing a direct proof of the validity of these conclusions. The open questions that require 
conclusive answers are: 

• Composition. Is the composition indeed dominated by protons, or is there a transition back to heavier nuclei at 
the highest energies? What is the cross section for pp interaction at high, > 100 TeV, energy? 

• Galactic- XG transition. At what energy does the flux become dominated by XG sources? 
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FIG. 10: GRBs and AGN are believed to be powered by black holes. The accretion of mass onto the black hole, through 
an accretion disk, releases large amounts of gravitational energy. If the black hole is rotating rapidly, another energy source 
becomes available: The rotational energy may be released by slowing the black hole down through interaction with the accretion 
disk. The energy released drives a jet-like relativistic outflow. The observed radiation is produced as part of the energy carried 
by the jets is converted, at large distance from the central black hole, to electromagnetic radiation. 

• Sources. Are the sources indeed transient? If so, are the sources GRBs or other transients? 

• Acceleration. Are UHECRs accelerated, as suspected, in collisionless (relativistic) shocks? A theory of such 
shocks based on basic principles is not yet available (e.g. [98| and refernces therein). 

In addition to the open questions listed above, the physics of the candidate UHECR sources is also not well 
understood. As we have shown, UHECR sources are required to produce very large power and are likely to be driving 
relativistic outflows, see Eqs. (JT]) and ©. These requirements suggest that the sources are powered by the accretion of 
mass onto black holes, as believed to be the case for GRBs and AGN. GRBs are most likely powered by the accretion 
of a fraction of a Solar mass on a ~ 1 s time scale onto a newly born Solar mass black hole El HI, I1H . Recent 
observations strongly suggest that the formation of the black hole is associated with the collapse of the core of a very 
massive star. AGN are believed to be powered by accretion of mass at a rate of ~ 1 Solar mass per year onto massive, 
million to billion Solar mass, black holes residing at the centers of distant galaxies 64]. As illustrated in figure ITUI 
the gravitational energy released by the accretion of mass onto the black hole is assumed in both cases to drive a 
relativistic jet, which travels at nearly the speed of light and produces the observed radiation at a large distance 
away from the central black hole. The models describing the physics responsible for powering these objects, though 
successful in explaining most observations, are largely phenomenological: The mechanism by which the gravitational 
energy release is harnessed to drive jets, the mechanism of jet collimation and acceleration, and the process of particle 
acceleration (and radiation generation), are not understood from basic principles. In particular, the answer to the 
question of whether the jet energy outflow is predominantly electromagnetic or kinetic, which has major implications 
to our understanding of the mechanism by which the jets are formed, is not known despite many years of photon 
observations. 

These open questions are unlikely to be answered by UHECR observatories alone. For example, given the un- 
certainties in the high energy pp cross section, it is not clear that studying shower properties would determine the 
primary composition. The composition could be constrained by an energy dependent anisotropy study (see § III CI) . 
However, the conclusions of such an analysis would depend on some assumptions regarding the sources [69j. In ad- 
dition, UHECR observatories are unlikely to identify the sources. Although they may provide a conclusive evidence 
for the correlation between the distribution of UHECR sources and that of matter in the local universe, and possibly 
discriminate between steady and transient sources (which may require large exposure that can be provided only by 
space-born detectors, see § III Gl) . this would not determine which type of objects the sources are. It should be em- 
phasized that electromagnetic observations are equally unlikely to resolve the open questions: despite many years of 
observations we are still lacking direct evidence for acceleration of nuclei in any astrophysical object, and fundamental 
questions related to the physics of the sources (e.g. the content of relativistic jets) remain unanswered. 

Thus, resolving the UHECR puzzles would require a "multi-messenger" approach, combining data from UHECR, 
7-ray and neutrino detectors. Neutrino astronomy is likely to play an important role in this context: detection of 
GZK neutrinos (see § IIIIB[) . combined with accurate measurements of the UHECR flux and spectrum, may allow us 
to determine the UHECR composition (and constrain the UHE pp and neutrino interaction cross sections); detection 
of high energy neutrino emission from electromagnetically identified sources may allow us to identify the UHECR 
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sources; Neutrino observations will provide new constraints on the physics driving the sources, which can not be 
obtained using electromagnetic observations, since they can escape from regions which are opaque to electromagnetic 
radiation (see § IIII CI for examples related to GRBs). 

Finally, it should be realized that if the UHECR sources are steady, identifying the sources by directly detecting 
their neutrino emission is highly improbable, due to the fact that the effective area of a 1 km 2 neutrino detector is 

3 x 10~ 4 km 2 at 10 3 TeV, ~ 10 -7 of the area of > 10 19 eV CR detectors (hence, neutrinos will not be detected unless 
the neutrino luminosity of the sources exceeds their UHECR luminosity by a factor of 10 3 ). In this case, identifying 
the sources will require a theoretical analysis combining electromagnetic, CR and neutrino data. 
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